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Abstract. The surface magnetic field configuration of the Ap star HD 118022 (78 Vir) has been reconstructed 
in the framework of the magnetic charge distribution (MCD) method from the analysis of Stokes IQUV spectra 
obtained using the MuSiCoS spectropolarimeter at Pic du Midi Observatory. Magnetically-sensitive Fe II lines 
were primarily employed in the analysis, supposing that iron is evenly distributed over the stellar surface. We 
show that the Stokes IQUV profile shapes and variations of 78 Vir can be approximately fit assuming a global 
magnetic field configuration described by a slightly decentered, inclined magnetic dipole of polar surface intensity 
approximately 3.3 kG. The derived inclinations of the stellar rotational axis to the line of sight i = 24 ±5° as well 
as to the magnetic dipole axis /3 — 124 ± 5° are in good agreement with previous estimations by other authors, 
whereas the sky-projected position angle* of the stellar rotation axis Q ~ 110° is reported here for the first time. 
In addition, several lines of Cr II and Ti II were studied, yielding evidence for non-uniform surface distributions 
of these elements, and magnetic field results similar to those derived from Fe. 
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1. Introduction 

78 Virginis (HD 118022) is a bright Ap star, and the first 
star other than the sun in which magnetic field was dis- 
covered (Babcock ll947fl . The longitudinal magnetic field 
of 78 Vir, as first observed by Babcock, is variable, and 
constantly negative. 78 Vir is also marginally variable 
in broadband light, in radial velocity, and shows varia- 
tions of its line profiles (Preston ©69). As suggested by 
Preston (196HJ), the observed properties of 78 Vir can be 
explained in the framework of the oblique rotator model 
(Stibbs 1950), in which the star rotates with a period of 
approximately 3.7 days and has a rotation axis forming a 
small angle with respect to the line of sight. 

The first comprehensive modelling of the surface mag- 
netic field structure of 78 Vir was performed by Borra 
(1980). Borra obtained high-resolution measurements of 
circular polarisation across the Fell A4520.2A line, which 
he attempted to reproduce using nine different magnetic 

Send offprint requests to: V. Khalack 
e-mail: khalakv@umoncton.ca 

* Q increases clockwise from the axis to the North Celestial 
Pole and relates to the azimuth angle O specified by Landolfi 
at al. ( T5331 as n = 360° - 6. 



field configuration models. As argued by Borra ©80), 
the profiles were well reproduced with an inclined dipole 
geometry that contains a moderate quadrupolar compo- 
nent. Nevertheless, he also pointed out that a decentered 
(a=0.2) dipole model provides essentially the same fit 
quality. Borra also concluded that the 78 Vir presents to 
an observer a remarkably uniform magnetic field over most 
of its visible disk, at most phases. It is unfortunate that 
due to the small inclination of the rotational axis (about 
25°) a significant part of the stellar surface remains hidden 
- a part which might have a less uniform field geometry. 

The analysis of broadband linear polarization (BBLP) 
measurements of several Ap stars (including 78 Vir) by 
Leroy et al. H199ri|) suggests that the magnetic fields 
of many Ap stars exhibit departures from the standard 
oblique rotator model assuming a pure dipole field ge- 
ometry. According to those authors, significant discrepan- 
cies between the BBLP observations and the "canonical 
model" results (Landolfi et al. I1993|) required the assump- 
tion of local departures from a dipolar field. In partic- 
ular, Leroy et al. l|199rj[) showed that the BBLP varia- 
tions of most stars could be reproduced assuming a dipole 
magnetic field geometry with slightly expanded field lines 
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over some parts of the magnetic equator. For 78 Vir, the 
observed BBLP variations do not show especially strong- 
departures from the dipolar case, although according to 
Leroy et al. ( 1 996 i a modified dipolar model provides a 
somewhat better fit to the BBLP curves for this star. 

In order to better constrain the magnetic field con- 
figuration of 78 Vir, in this paper we undertake a more 
detailed modelling of the surface magnetic field structure 
based on the analysis of high resolution line profiles in all 
4 Stokes parameters. In Sec. we discuss the properties 
of obtained Stokes IQUV spectra, and in Sec.|3]we sum- 
marise the fundamental characteristics of 78 Vir. In Sec. 01 
we describe the features of the magnetic field modelling 
framework, and derive the global magnetic field configura- 
tion of the star, comparing observed and computed Stokes 
profiles for various spectral features. In Sec. |3]we discuss 
the derived global magnetic field characteristics and their 
agreement with earlier studies. 

2. Observations 

Spectropolarimetric observations of 78 Vir were obtained 
in 1997 February, 1998 February and 1999 January using 
the 2m Telescope Bernard Lyot at Observatoire du Pic du 
Midi (Wade et al. I2000a|) . The MuSiCoS cross-dispersed 
echelle spectrograph (Baudrand & Bohm 1992) and ded- 
icated polarimeter module (Donati et al. 1999) were em- 
ployed for the observations. 

The MuSiCoS spectrograph is a table-top instrument, 
which allows the acquisition of a stellar spectrum in a 
given polarization state (Stokes V, Q or U) throughout the 
spectral range from 4500 to 6600 A with a resolving power 
of about 35 000, in a single exposure. The spectrograph is 
fed by a double optical fibre directly from the Cassegrain- 
mounted polarimeter. The optical characteristics of the 
polarization analyser, as well as the spectropolarimeter 
observing procedures, are described in detail by Donati 
et al. l|1999f) . Observing details specific to the acquisition, 
reduction and analysis of the 78 Vir spectra are provided 
by Wade et al. l|2000a|l . 

The journal of spectropolarimetric observations is re- 
ported by Wade et al. I|200()a|l . The 52 Stokes V, Q and U 
spectra, obtained on 18 different nights, cover the whole 
rotational period of 78 Vir approximately uniformly, and 
provide an average S/N of about 370. 

3. Fundamental parameters of 78 Vir 

78 Vir was classified by Cowley et al. (jT069 ) as 
AlpCrSrEu. Adelman l)1973al I1973b|l performed both a 
line identification and abundance analysis of this star. 
The distance to 78 Vir d = 48 ± 15 pes and its radius 
R = 1.77 ± O.68i? are derived by Monier iTH^ using 
the Infrared Flux Method. On the other hand, taking into 
account the mean angular diameter 9 a = 0.343 milliarcsec 
obtained by Monier l|1992[l for this star and its distance 
d = 56.2 ± 2.5 pc derived from the Hipparcos parallax 



Table 1. Fundamental parameters of 78 Vir. 



Parameter 


Value 


Reference 


Spectral type 


AlpCrSrEu 


Cowley et al. (I1969t 


Age 


3.0±°{ 7 3 x 10 8 y 


this work 


Distance 


56.2±2.5 pes 


Hipparcos (ESAll99?ll 


Period 


3. d 7220±0. d 003 


Preston fTMfll 




9200±290K 


Monier l|1992(l 


log5 


4.50± 0.25 


ibidem 


M v 


1.18±0.10 


this work 


L* 


27.3±2.5 L Q 


this work 


R+ 


2.O6±O.17J? 


this work 


M+ 


2.18±O.O6M 


this work 


V c sin i 


12±1 km s" Y 


this work 


i 


25°±5° 


Leroy et al. iTf^pl 


P 


120°±5° 


ibidem 



fESA lTMTIl . we find a somewhat larger value for the stel- 
lar radius R = 2.06 ±0.17 R@. Using the Hipparcos vi- 
sual magnitude and parallax we have found the absolute 
visual magnitude of 78 Vir, M v = 1.18 ± 0.10. Taking 
into account the bolometric correction (Flower ^996), 
which corresponds to Monier's (1992) 9200 ± 290 K ef- 
fective temperature and the bolometric zeropoint correc- 
tion BCv=-0.07, which is obtained assuming M b ® =4.74 
for the Sun (Bessell et al. 1998|, we can find the ab- 
solute luminosity for 78 Vir L* = (27.3 ± 2.5)L Q . For 
78 Vir Monier 1|1992J) also estimates the integrated flux 
/* = (2.7±0.27) x 10~ 7 erg sm^s^ 1 . Taking into account 
the distance to the star, this provides the luminosity L* = 
(26.5 ± 2.8)Lq. Finally, employing the stellar radius and 
the effective temperature (see Table we obtain a third 
estimate of the absolute luminosity ~ 27 .4 ± 4.5Lq. 
All of these values are in good mutual agreement. The 
luminosity = (27.3 ± 2.5) Lq together with the effec- 
tive temperature (see Tableful allow us to find the stellar 
mass M* = 2.18 ± O.O6M and age 3.0±5'.3 x 10 8 years 
(see Fig. by spline interpolation (Sandwell 1987) in 
the model evolutionary tracks of Schaller et al. (1992) for 
metallicity Z=0.02. The derived age suggests that 78 Vir 
has completed approximately 37% of its main sequence 
life. The mass and radius provide logg = 4.16 ± 0.07. 
This value is marginally inconsistent with that of the 
Monier (1992), but it is in good agreement with the value 
log g — 4.20 reported by King et al. (£2003). A recent in- 
vestigation of the Ursa Major stream characteristics per- 
formed by King et al. (|2TjTj3) derived Z = 0.016 - 0.02 and 
age (5 ± 1) x 10 8 years for the stream. They found that 78 
Vir is a " certain member" of the stream from photometric 
data, while the kinematic characteristics of the star argue 
for its "probable non-membership". The derived age for 
78 Vir is marginally inconsistent with the UMa stream 
age. The published metallicity of 78 Vir is log(Fe/H")* — 
log(Fe/H) Q =+°:H (Cayrel de Strobel et al. 1197)71) and 
it is not clear if this star truly has a solar metallicity. For 
this reason the derived uncertainties of the mass and age 
may be somewhat larger than indicated here. 
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Babcock (1947) observed that the longitudinal mag- 
netic field of 78 Vir was variable, but always negative. 
Analysing Babcock's data along with his own measure- 
ments, Preston l|1969jl determined the periodic character 
of the magnetic field variability and derived the ephemeris: 

JD (magnetic maximum) = 2434816.9 + 3. d 7220 • E. (1) 

Several authors have confirmed this period on the basis of 
further magnetic field measurements (Wolff & Wolff ESI 
Wolff & Bonsack[M2 Wolff [M| BorraQjEOI Borra fc 
Landstreet HMQl Landstreet 113521 Wade et al. 12000 b|l . 
while Landstreet (private communication) has estimated 
the accuracy of the period determination to be ±0. d 003. 
Preston 119691) also found that the crossover effect is 
very pronounced and is present throughout most of the 
magnetic cycle. This feature is particularly remarkable 
at phase 0.85 (Wolff & Bonsack 113721 see also Wade et 
al. l2000al) . 

Applying the oblique rotator model with the dis- 
torted dipole approximation, Leroy (private communica- 
tion) found that they must adopt a period 3.7218 days 
in order to fit the observed BBLP variations. As dis- 
cussed in Sect. ^ their model, using local modifications 
of the axisymmetric magnetic field imposed near the 
magnetic equator, differs only mildly from a dipole con- 
figuration. The derived rotational axis inclination and 
dipole obliquity i = 25° and (3 — 120° provide a good 
match to the longitudinal field and linear polarisation 
variations. Nevertheless, the period inferred is somewhat 
shorter than that obtained by other authors (Eq. , and 
Wade et al. (2000b) reported that this period is only 
barely consistent with the available longitudinal field mea- 
surements. 

78 Vir also shows photometric variability in the visible 
and the infrared with the same period (Eq.f^). Combining 
the visual flux variability in the uvby bands, obtained by 
Wolff & Wolff l)1971|l . with their own new observations, 
Catalano & Leone l|1994[l have tried to improve the period 
determination. Their analysis results in the new ephemeris 

JD (y min.) = 2434816 + (3. rf 722084 ± d 000042)E, (2) 

which was recently confirmed by Leone & 
Catanzaro l|2001|l . These authors have shown with 
the help of Hipparcos photometry that the light vari- 
ations of 78 Vir are not purely sinusoidal. Taking into 
account all archival magnetic field observations as well 
as their own data, Leone & Catanzaro (2001) have 
shown that the longitudinal field measurements are in 
phase only when the 3.722084 day period is adopted. 
Moreover, by adopting the 3.7218 day period suggested 
by Leroy (private communication), they find a 0.12-cycle 
phase shift between the Hipparcos light curves and the 
other light curves. 

Monier (1992J has derived for 78 Vir from the sim- 
ulation of the energy distribution in the spectral range 
from 1200A to 22000A an effective temperature T c g = 



1.8 



1.6 



1.2 



ZAMS\ 










\ 2.5/ 










2.18/J, 

\ 










1.7^--— " 



4.1 4.05 4 3.95 3.9 3.85 3.8 



logT e „[K] 

Fig. 1. The position of 78 Vir on the H-R diagram 
is compared with the theoretical evolutionary tracks of 
Schaller et al. I|1992|l . computed for solar metallicity and 
masses M = 1.7, 2.0 and 2.5 Mq. The age 3.0 x 10 8 years 
and the mass M* = 2.18 Mq of the star are estimated by 
interpolation of the theoretical evolutionary tracks. The 
dashed line shows the isochrone, which corresponds to the 
age of 78 Vir (3.0 x 10 8 y). 

9200 ± 290 K, a surface gravity logg = 4.5 and a pho- 
tospheric metallicity [M / H] = 10. In addition, using a 
model atmosphere with these values of log g and [M/ H] , 
he has found that the effective temperature T c g = 9300K 
provides the best description of the spectral energy dis- 
tribution in the region from 1200A to 8000A at rota- 
tional phase 0.0, while T c ff = 9200K seems to be the best 
effective temperature, at phase 0.5. These temperatures 
are somewhat lower than most previous estimates, rang- 
ing from 9700K to 10700K (Mihalas & Henshaw ITUCT 
Wolff HM7I Jugaku & Sargent Il968f) . probably because 
the models used by these authors were unblanketed and 
calculated for solar abundances. 

The fundamental parameters of 78 Vir are summarised 
in Table □ 

4. Modelling of the Stokes IQUV spectra 

According to the discussion of Sect. [31 an ATLAS9 model 
fKurucz !T994|l with parameters T c g = 9250K, \ogg = 4.5, 
[M/H]=0 and microturbulent velocity v t =0 km s -1 suc- 
cessfully approximates the stellar atmosphere of 78 Vir. 

We describe the structure of the surface magnetic field 
of 78 Vir (assumed to be a rigidly rotating and spherically 
symmetric star) using the magnetic charge distribution 
(MCD) method (Gerth et al. lTMTI Khalack et al. EDOIl). 
Originally, the MCD method considered a system of spa- 
tially separated point field sources with virtual magnetic 
charges in the stellar interior. In order to provide zero 
magnetic flux through the stellar surface the sum of "mag- 
netic charges" should be kept to zero. These sources 
produce a magnetic field whose potential at each point 
on the stellar surface is specified by the superposition 
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Table 2. Here the individual columns specify the ion, 
the wavelength, the quantum number J and the Lande 
factor for lower and upper atomic levels, the observed 
Stokes / profile depth and the line scaling factor (LSF) 
(see Sect. I4.3[) . derived from the comparison of simulated 
equivalent width of the Stokes Q and U Zeeman profiles 
with the BBLP data (Leroy [1995). Asterisks marks the 
LSF calculated from Stokes Q and U profiles simulated 
for a field structure that is derived from only Stokes / and 
V profile variability (see Table. |3J). 



Ion 


A, A 






J up 


ffup 


-^obs 


LSF 


Fen 


4620.52 


3.5 


1.21 


3.5 


1.40 


0.33 


0.035 


Fen 


4635.32 


2.5 


1.20 


3.5 


1.13 


0.32 


0.036* 


Fei 


4635.85 


1.0 


1.49 


2.0 


1.89 






Fen 


4923.93 


2.5 


2.00 


1.5 


2.40 


0.58 


0.079 


Til 


5017.95 


4.0 


1.50 


5.0 


1.18 






Cri 


5018.15 


2.0 


1.16 


3.0 


1.09 






Fen 


5018.44 


2.5 


2.00 


2.5 


1.87 


0.60 


0.102 


Cm 


5018.84 


3.5 


1.24 


3.5 


1.42 






Fen 


5100.61 


4.5 


1.54 


4.5 


1.34 






Fen 


5100.66 


4.5 


1.31 


3.5 


1.40 






Fen 


5100.73 


4.5 


1.54 


5.5 


1.36 






Fen 


5100.85 


1.5 


0.80 


1.5 


0.72 


0.47 


0.083* 


Fei 


5168.90 


3.0 


1.50 


3.0 


1.75 






Fen 


5169.03 


2.5 


2.00 


3.5 


1.70 


0.59 


0.178 


Fei 


5169.30 


4.0 


1.26 


3.0 


1.32 






Fen 


5197.48 


2.5 


1.20 


1.5 


0.72 






Fen 


5197.58 


2.5 


0.57 


1.5 


0.44 


0.39 


0.073 


Fen 


5362.74 


3.5 


1.15 


4.5 


1.24 






Fen 


5362.87 


4.5 


1.15 


3.5 


1.40 


0.45 


0.112* 


Cri 


5362.96 


3.0 


1.33 


3.0 


1.06 






Fen 


5362.97 


3.5 


1.44 


4.5 


1.34 






Crn 


5363.88 


3.5 


1.71 


3.5 


1.39 






Fen 


6247.35 


1.5 


0.40 


2.5 


0.62 






Fen 


6247.56 


2.5 


1.33 


1.5 


1.72 


0.36 


0.036* 


Fen 


6432.68 


2.5 


2.00 


2.5 


1.65 


0.24 


0.026 


Fen 


6516.08 


2.5 


2.00 


3.5 


1.58 


0.25 


0.032 



of potentials of individual sources (Gerth et al. [1997 
Gerth & Glagolevskij I3UUT1 12TM|> . Since the number of 
sources is usually more than one, we actually operate with 
a system of several magnetic dipoles located in the stellar 
interior. In order to minimize the number of free model pa- 
rameters, the most convenient way is to consider a system 
of two sources, which mathematically formulate a mag- 
netic dipole (Khalack et al. J2001). When the dipole is 
centered on the stellar centre and the dipole size is much 
smaller than the stellar radius fKhalack l2()02|) . the MCD 
model transforms to the conventional model of a symmet- 
ric inclined magnetic rotator fStibbs lT95()(l . Otherwise, we 
deal with a more complex magnetic field configuration. The 
mathematical verification of the MCD model as well as 
the procedure of specification of the angle (3 between the 
rotational and magnetic dipole axes, the coordinates and 
the field strengths of the positive and negative magnetic 
poles on the basis of the derived free model parameters is 
described in detail by Khalack et al. (p!003 ) . 



4.1. Procedure 

From previous analyses of the Stokes IQUV spectra of 
78 Vir, it has been found by Wade et al. (|2()()()a[> that 
the three strong Fen lines A4923.93A, A5018.44A and 
A5169.03A show the strongest Zeeman signatures of any 
lines in the optical spectra of CP stars. The real mag- 
netic field structure of 78 Vir is expected to differ only 
marginally from a centered magnetic dipole (based on pre- 
vious modeling efforts: Borra 1980). The most attention is 
therefore paid to the variability of the Stokes / and V pro- 
files, which contain the most information about the global 
surface magnetic field configuration. The Stokes I varia- 
tion provides information about the surface distribution 
of the chemical abundance, as well as the stellar radial 
velocity V t and projected rotational velocity L c sinL The 
Stokes V profiles are sensitive to the global field configura- 
tion, in particular the lower-order multipolar components. 
The Stokes Q and U profiles also provide some constraint 
on the global field morphology, but are most sensitive to 
the smaller-scale structure of the field. Unfortunately, the 
available Stokes Q and U profiles have low relative S/N 
ratio (Wade et al. l2(JU(Ja|) . and in this study they are used 
primarily to check the results of the Stokes I and V profile 
simulations and to determine the sky-projected position 
angle of the rotational axis, ft. According to Khalack et 
al. <|2()()1II2()()3|I this angle (0 < ft < 360°) is counted clock- 
wise from the rotation axis to the North Celestial Pole in 
the plane of the sky, and is related to the azimuth angle 
6 specified by Landolfi et al. l|M)3) by ft = 360° - 9. 
The opposite counting of the position angle in the MCD 
model is compensated for by the negative sign in definition 
of the i? y -component of the local magnetic field. Hence a 
right-handed reference frame is applicable. 

We have examined the spectra for additional lines with 
prominent polarisation signatures, in order to compile a 
list of lines that are especially appropriate for this task. 
To perform the line identification we use the VALD-2 re- 
sources (Kupka et al.:1999 Ryabchikova et al. H999|l . The 
adopted list contains Fen, Fei, Cm, Cri, Tin, Til and 
Mgn lines, but for the present study we will concentrate 
primarily on the strongest Fell lines. The main reason 
for this is that Fe is presumably almost uniformly dis- 
tributed over the stellar surface (see Sect. 14. 2|) and in this 
way we exclude from the fitting procedure a model of the 
surface abundance distribution. Table [21 presents the final 
list of Fe II lines used in the simulation procedure, together 
with lines of some other elements that are responsible for 
blends. 

Some of the Fell lines listed in Table have no 
detectable features in the Stokes Q and U spectra. 
Nevertheless, they are comparatively strong lines (in the 
Stokes / spectra) with clear variability in the Stokes V 
spectra, and are analysed without the linear polarization 
data in order to check our final results. 

All line profiles are simulated using the Zeeman 2 
polarised spectrum synthesis code (Landstreet [1988 
Wade et al. I2001fl . The code has been modified to in- 
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elude a magnetic field described within the framework of 
the MCD method (Khalack et al.[2003), and to allow for 
an automatic minimization of the model parameters us- 
ing the downhill simplex method (Press et al. 11992(1 . The 
relatively poor efficiency of the downhill simplex method, 
requiring a large number of function evaluations, is a well- 
known problem. Repeating the minimization routine 3-^4 
times in the vicinity of a supposed minimum in the pa- 
rameter space allows us to check if the method converges 
to a global minimum. In our case, this technique requires 
a comparatively long computational time due to the large 
amount of analysed observational data and the large num- 
ber of free parameters. 

The free model parameters employed in the line pro- 
file simulation can be divided into two groups: those that 
describe the magnetic field structure, and those that de- 
scribe the stellar geometry and atmospheric characteris- 
tics. The first group includes the parameters Q T , a\ and 
CL2 (which specify the modulus of "magnetic charges" and 
their distance from the center of the star, expressed in 
the units of stellar radius), and the parameters \%,5i 
and A2 , 62 (which specify the spherical coordinates of the 
"charges" in the rotational reference frame of the star). 
Meanwhile, the second group includes the parameters fl 
(the position angle) , i (the rotational axis inclination with 
respect to the line of sight), V r (the heliocentric radial 
velocity), V c sini (the projected rotational velocity), and 
\og(N x /N to t) (the abundance(s) of the element(s) forming 
the line to be modeled). Instead of V e smi we can obviously 
use only V e as a free parameter, but V e sin i is preferable 
here, because it allows us to compare our estimate of this 
parameter directly with the results of other authors. The 
minimum number of free model parameters is 11 for the 
MCD model if we do not take into account the linear po- 
larization data, and 12 if we do. Sometimes, to obtain a 
good fit to the line profiles we need to include in the sim- 
ulation lines of other chemical elements (see Table 0). In 
this case the number of free model parameters grows with 
the number of chemical elements under investigation, but 
usually does not exceed 15 parameters. 

Initially we specify arbitrarily the values of the free 
model parameters in the range of their possible values and 
simulate the Stokes IQUV spectra at the phases of the 
observations. We then compare the simulated profiles with 
the observed spectra, and calculate the reduced x 2 > which 
we adopt as a measure of the fit quality. The expression for 
the x 2 function, which reflects the agreement between (for 
example) the simulated Stokes / (1^ {(fii)) and observed 
Stokes / {I° b \.) spectra is given by: 



x 2 




(3) 



where cr[/°^ s ] corresponds to the measurement errors, (fii 
is the rotational phase of the observation, Nj specifies the 
number of spectra, while Ni is the number of pixels in 
each analysed line profile. In fact, although the simulated 



spectra are calculated with approximately the same re- 
solving power as the observed spectra, this does not pro- 
vide a direct coincidence of wavelengths in the simulated 
and observed spectra. Therefore, during the x 2 function 
evaluation the simulated spectral intensity at the exact 
observed wavelength is calculated using a linear interpo- 
lation. 

The scale of variability and the measurement errors 
are different for each of the four Stokes IQUV spectra 
(Wade et al. 2000a). In order to balance (from a statisti- 
cal point of view) the information flow from each of the 
four Stokes spectra, the respective \ 2 contributions are 
weighted. 

The weights are derived from the comparison of the 
best fit values of Xy 1 Xq 1 Xu with X\ (when the ordinary 
X 2 function is minimized) and provide approximately the 
same weighted value of these functions for the majority of 
analysed lines (see Table The weighted xt functions is 
specified in the following way: 



1 



[x? 



(4) 



(5) 



if we work only with the Stokes / and V spectra. The value 
of the Xw function is reduced throughout the minimization 
procedure, varying the values of all free model parameters. 
When it reaches its global minimum, we can also calculate 
the ordinary x 2 function: 



1 



X=jbCi+Xv + XQ+Xu]- 



(6) 



We begin by simulating the strong Fen A4923.927 A, 
A5018.44 and A5169.03 lines, for which the minimization 
process has been repeated for 13-^15 times starting from 
different locations in the free model parameter space in 
order to avoid local minima. In the case of 78 Vir, the 
downhill simplex method converges to four global min- 
ima, which are caused by the decentered dipolc model 
symmetry in the MCD method. Two global minima cor- 
respond to the parameter sets i,(3,fl and i,(3, 180° + fl 
which differ only by the angle f2 and provide exactly the 
same value of the x^,-function. Models with the oppo- 
site direction of stellar rotation (180° — i, 180° — j3, or 
180° - i, 180° - (3, 180° + n) do not result in minima. The 
other configurations (i, 180° -j3, n and i, 180° -0, 180°+f7) 
correspond to models in which the positive magnetic pole 
faces the observer. These configurations are not applica- 
ble to the case of 78 Vir, where we always see the neg- 
ative magnetic pole, while the positive pole is partially 
visible just for phases close to ip = 0.0 (Babcock 119471 
Borra H"9"8Ti|l . The other two global minima have parameter 
sets which differ from the previous ones by the location of 
the "magnetic charges" . The two "magnetic charges" can 
be located in two different stellar hemispheres (separated 
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Table 3. Results of Fen lines simulation for T e s = 9250K, log g = 4.5 and v t —Q km s^ 1 . The first column indicates 
the free parameters, while the other columns specify the parameter values for the given line profile. The last column 
provides the averaged estimation errors for each parameter. The first 6 rows show the fit quality for the each analysed 
Stokes profile (Eq.[3J), weighted and unweighted (Eqs.0EJ| X 2 ~f unc ti° n - The following 12 rows show the best fit values 
of the free model parameters, while the final 9 rows provide the characteristics of the magnetic dipole poles, which are 
derived from the free parameters. 



Line, A 


4620 


4635 


4923 


5018 


5100 


5169 


5197 


5362 


6247 


6432 


6516 




Xi 


2.92 


10.91 


8.96 


6.97 


7.40 


13.10 


16.06 


5.54 


2.27 


8.83 


5.94 




4Xv 


6.14 


7.86 


9.14 


8.69 


6.97 


9.09 


8.71 


6.64 


5.00 


12.24 


7.39 




6Xq 


7.75 




7.73 


5.76 




8.23 


6.35 






12.15 


9.30 




7.42 




10.04 


6.82 




10.79 


6.80 






9.99 


6.19 






6.06 


9.39 


8.97 


7.06 


7.19 


10.30 


9.45 


6.09 


3.64 


10.81 


7.20 




x 2 


1.75 


6.44 


3.55 


2.81 


4.57 


4.64 


5.11 


3.60 


1.76 


3.90 


2.59 




Qr, kG 


180 


254 


303 


202 
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221 


151 


192 


398 


173 


289 


50 




8.4 


6.7 


4.3 


6.7 


5.3 


7.3 


12.6 


9.7 


5.4 


8.9 


6.4 


1.2 


Ai 


13° 


24° 


34° 


41° 


33° 


15° 


-4° 


19° 


69° 


23° 


20° 


5° 


Si 


-45° 


-46° 


-52° 


-51° 


-60° 


-30° 


-40° 


-44° 


-39° 


-45° 


-39° 


8° 


a 2 ,l(T 3 
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5.2 


4.3 


2.5 


1.6 


3.7 


5.1 


3.8 


2.7 


0.3 


A 2 


118° 


125° 


144° 


135° 


145° 


85° 


111° 


101° 


113° 


96° 


111° 


10° 


52 


-16° 


-9° 


-9° 


-14° 


-12° 


-11° 


-11° 


-26° 


-11° 


-27° 


-11° 


10° 





109° 




109° 


109° 




109° 


102° 






164° 


128° 


17° 


i 


22° 


24° 


22° 


23° 


25° 


21° 


29° 


27° 


26° 


29° 


24° 


5° 


logGFe/iVtot) 


-3.36 


-2.82 


-3.26 


-3.24 


-3.20 


-3.12 


-3.25 


-2.86 


-3.17 


-2.97 


-3.50 


0.22 




-8.78 


-8.42 


-7.82 


-7.51 


-8.69 


-7.23 


-7.93 


-8.25 
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-6.57 
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1.18 


V e sin i 


10.9 


11.1 


11.9 


12.0 


11.5 


12.6 


11.9 


13.0 


11.9 


13.7 


11.5 


0.9 


(3 


125° 


124° 


122° 


120° 


123° 


118° 


129° 
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123° 
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1.6 


a,l(T 3 


4.5 


3.8 


2.6 


0.7 


3.4 


3.4 


6.2 


4.4 


2.1 


3.8 


3.4 


1.6 


B v , kG 


3.25 


3.95 


3.47 


3.18 


3.13 


3.02 


4.00 


3.45 


3.42 


2.69 


3.92 


0.7 


\p 


-9° 


-10° 


-6° 


-7° 


-10° 


-8° 


-11° 


-8° 


-11° 


-8° 


-8° 


5° 


s p 


-35° 


-33° 


-30° 


-30° 


-33° 


-28° 


-37° 


-36° 


-34° 


-36° 


-32° 


12° 


B n , kG 


-3.18 


-3.82 


-3.46 


-3.16 


-3.12 


-2.95 


-3.85 


-3.36 


-3.41 


-2.62 


-3.86 


0.7 


X n 


170° 


170° 


174° 


173° 


170° 


172° 


168° 


171° 


169° 


171° 


172° 


5° 


5n 


35° 


33° 


30° 


30° 


33° 


28° 


37° 


35° 


33° 


36° 


32° 


12° 



by the equatorial plane), or in a single hemisphere. If they 
are not significantly shifted from the stellar centre, they 
can form the same dipole axis and will produce almost the 
same surface magnetic field configuration. The differences 
between the two models can be revealed by the sensitivity 
of Xw function on the parameter values, but the minima 
will still exist. Preference is given to the deepest global 
minimum, obtained for the model with the two "magnetic 
charges" located in the same stellar hemisphere. For the 
other lines the minimization process was performed only 
3-^4 times (although if we tested the contribution of blends 
to the analysed profiles, we ran the minimization routine 
several more times). Supposing that all analysed lines con- 
tain the signatures of the same magnetic field structure, 
we chose the initial locations in parameter space not far 
from the parameter set obtained from the strong line anal- 
ysis. 

In order to evaluate the fit errors (and therefore the 
uncertainties on the derived free parameters) , we calculate 
deviations of the simulated profiles produced as a result 
of small variations of each of the free parameters, thus 
introducing a small shift along one axis in the x 2 hyper- 
space from the point of the function minimum value. Using 



this procedure, and taking into account the uncertainties 
of the observational data and the obtained minimal value 
of x 2 -function, we can estimate the errors of the best-fit 
parameters. 

4.2. Results 

Our initial assumption that iron is uniformly distributed 
over the surface of 78 Vir appears to be valid, given that 
the analysed Fe n lines reveal no significant variability of 
the Stokes / profiles with rotational phase. Each Fen line 
(or group of lines) shown in Table ^ was analyzed inde- 
pendently of the others, using a stellar atmosphere model 
with T cff = 9250K, logy = 4.5, v t =0 km s _1 . For each 
best-fit simulation the derived free model parameters are 
given in Tabled During the simulation of some Fen lines 
the Stokes Q and U profiles were not taken into account, 
and consequently the angle f2 is not determined for those 
lines. 

The other data in Table 01 are derived from the free 
model parameters (Khalack et al. 2003). Here (3 defines 
the angle between the magnetic dipole axis and the stel- 
lar rotation axis (the angle exists if these two axes cross 
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a) 7/7 c 



b) V/h 



c) Q/Ic 



d) U/Ic 




5018 



5019 



5018 



5019 



5018 



5019 



5018 



5019 



Fig. 2. Strong lines. Observed (solid lines) and simulated (dashed lines) Stokes I (a), V (b), Q (c) and U (d) line 
profiles for Fen A5018.44A line. The vertical bars show the size of the respective observational errors. The observed 
and simulated profiles are shifted by 0.15, 0.045, 0.02 and on 0.015 along the vertical axes for Stokes I, V, Q and U 
profiles respectively. 



each other). The distance of the magnetic dipole center 
from the center of the star and one-half of the magnetic 
dipole size are represented by variables clq and a, respec- 
tively. The variables (B p , X p , S p ) and (B n ,X n ,6 n ) specify 
the location of the positive and negative magnetic poles 
at the stellar surface and the respective strength of the 
magnetic field. 

4.2.1. The strong lines 

The Fen A4923.927, A5018.44 and A5169.03 lines are the 
most prominent Fe lines in the spectrum of 78 Vir, and 
show the clearest Stokes Q and U signatures. The first pro- 
file is formed essentially by the single line Fen A4923.927 
(see Table|2J) and appears to contain no important blends. 
The agreement between the observed and simulated data 
for Fen A4923.927 is similar to that shown in Fig. (for 
Fe II A5018). Independent Stokes I spectra (each corre- 
sponding to a slightly different phase) were obtained with 
each of the Stokes V, Q and U spectra, and although only 



one set of Stokes / profiles are presented in the figure, all 
were taken into account during the calculation of the Xi~ 
function (Eq.EJ. The best-fit values of the Stokes IVQU 
X 2 -functions for the Fell A4923.927 line are given in the 
third column of Table |21 

The second line profile is composed mainly of the Fe II 
A5018.44A line, but is also contaminated by contributions 
from the Til A5017.95, Cri A5018.15 and Cm A5018.84 
lines. It seems that the Cm A5018.84A is responsible for 
the blend in the red wing of the Stokes I and V profiles 
(Fig. [2J. Supposing a uniform Cr distribution, this blend 
is well fit for a Cr abundance of log Cr/iV tot =-3. 38 dex. 
The best fit x 2 values for the Stokes IV QU profiles of Fe n 
A5018.44 are given in the fourth column of Table |3| 

In the case of Fen A5169.033, the profile is blended 
by the weaker Fei A5168.898 and A5169.296 lines (see 
Table EJ. According to Kochukhov et al. H2004[> log#/=- 
0.786 given in the GRIFON list for the Fei A5169.296A 
line is too high to match the solar spectrum, and they 
recommend to use a decreased log g/=-2.15. This value 
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a) J/7 C 



b) V/h 



c) Q/h 



d) U/h 




Fig. 4. 



6516 6516.5 

The same as at the Fig. |2| 



6516 6516.5 

but for Fen A6516 



6516 6516.5 



6516 6516.5 



> 
< 




1.25 



Fig. 3. The difference in radial velocity between the 
observed and calculated Stokes I profiles for the Fell 
lines A4923.927A (crosses), A50f8.44A (open circles) and 
A5169.033A (open squares). The mean radial velocities 
V r =-7.82 km s _1 , -7.51 km s" 1 , -7.23 km s" 1 obtained 
from the simulation of Fen A4923.927A, A5018.44A, 
A5169.033A lines respectively (see Table 01, have been 
taken into account. 



is employed here for the respective profile simulation and 
provides much better agreement with the observed data. 
The behaviour of the best fit Stokes IVQU profiles with 



stellar rotational phase for this particular line is very sim- 
ilar to that shown in Fig. for Fen A5018.44. The sixth 
column of Table|3presents the best fit parameters and the 
model characteristics for Fen A5169.033. 

The concordance of the observed and simulated 
Stokes I and V profile variations as a function of stellar 
rotation is extremely good for these lines. On the other 
hand, the fit to the Stokes Q and U profiles is only ap- 
proximate. The simulated profiles show generally the same 
intensity, qualitative structure and variability as the ob- 
servations. However, even with these relatively noisy data, 
it is clear that the model does not reproduce the observa- 
tions within the errors at some phases (for example, see 
phases 0.1174, 0.5339 and 0.9851 in Fig. 2). 

In order to verify the assumption of a uniform surface 
distribution of iron in the atmosphere of 78 Vir, we have 
calculated the difference in radial velocity between the ob- 
served and simulated Stokes I profiles for the strong Fe n 
lines A4923.927, A5018.44 and A5169.033A (see Fig.|3J). 

The mean (averaged for all the available observational 
phases) radial velocities obtained from the simulation are 
taken into account in the calculation of the velocities of 
the respective simulated profiles. Fig.|3]shows that the de- 
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rived differences in V Y almost coincide for all three lines. 
A moderate disagreement may exist only for the data ob- 
tained in the vicinity of rotational phase </?=0.5, when the 
negative magnetic pole is most visible. The variation of 
radial velocity, which shows a reasonably coherent varia- 
tion with phase from about -1 to +1 km/s, may reflect a 
mildly non-uniform distribution of Fe, unmodelled struc- 
ture in the magnetic field, or other unaccounted-for phys- 
ical processes in the stellar atmosphere. 

4.2.2. The moderate strength lines 

From the list of the Fen lines selected for analysis 
(see Table the A4620.52, A5197.58, A6432.68 and 
A6516.08 lines are weaker than the A4923.927, A5018.44 
and A5169.03 lines. Due to the weaker polarised signal, a 
few spectra with comparatively large observational errors 
have been excluded from the simulation (phases 0.5757, 
0.5847 and 0.5914). 

The line Fe n A4620.521 is primarily responsible for the 
formation of the observed profile. No blends were taken 
into account during the simulation in this case. The Fe n 
A5197.577 line is blended by the weak Fen A5197.48 line, 
but provides the main contribution to the observed profile. 
The behaviour of the best fit Stokes IVQU profiles with 
stellar rotational phase for these lines is very similar to 
that shown in Fig. H for Fen A4620.52lA. The first and 
seventh columns of Table 01 present the best fit parameters 
and the model characteristics for Fen lines A4620.52lA 
and A5197.577A, respectively. 

The Fen A6432.68A line is also essentially unblended. 
The agreement between the observed and simulated 
Stokes IVQU data for this particular line is very simi- 
lar to that presented by Fig.H The Fen A6516.08A line 
is the primary contributor to the formation of the cor- 
responding observed Stokes IVQU profiles. The agree- 
ment of the best-fit simulation with the observed data 
for Fen A6516.08A line is shown at Fig. 01 The tenth and 
eleventh columns at the Tables 01 contain the best fit pa- 
rameters and model characteristics for the Fen A6432.68A 
and A6516.08A lines, respectively. The best fit of the sim- 
ulated Stokes IVQU profiles for the Fe n A6516.08A line is 
statistically better than that obtained for Fen A6432.68A. 

4.2.3. Separate analysis of the Stokes I and V spectra 

The other Fen lines A4635.316A, the A5100A- group, 
A5362.87A and A6247.56A show no variability in the 
Stokes Q and U spectra and hence only the Stokes / and 

V spectra are taken into account during the simulation. 
These lines are selected for the analysis in order to check 
the resulting magnetic field structure of 78 Vir on the basis 
of a more thorough line list. 

The Fen A4635.316 line is the main contributor to 
the formation of the corresponding observed Stokes / and 

V profiles, although it is blended by the comparatively 
weak Fei A4635.846A line. The agreement between the 



a) I/I c 



b) V/Ic 




5100.5 5101 



5100.5 5101 



Fig. 5. Observed (solid lines) and simulated (dashed lines) 
Stokes I (a) and V (b) profiles for Fen A5100.607A, 
A5100.664A, A5100.727A and A5100.852A lines. The ver- 
tical bars show the size of respective observational errors. 
The observed and simulated profiles are shifted on 0.15, 
and on 0.045 along the vertical axes for Stokes I and V 
profiles correspondingly. 



observed and simulated data for this particular line is sim- 
ilar to that shown in Fig. 01 The superposition of the Fen 
A5100.607, A5100.664, A5100.727 and A5100.852 lines is 
responsible for the formation of Stokes I and V profiles 
at ~ 5100.7 A (see Table 0). Fig. shows the simulated 
Stokes I and V profiles, which fit well the observed pro- 
files (see also the respective % 2 value in the fifth column 
of Table 0J). 

In the region of the Fen A5362.87 line the observed 
Stokes I and V profiles are formed mainly by this line 
and by a weak contribution from the Fen A5362.74, 
A5362.98, Cri A5362.87 and Cm A5363.88 lines. The Fen 
A6247.557 line also provides the main contribution to the 
observed Stokes / and V profiles and is blended by the 
Fen A6247.35 line. The best fit simulated data show al- 
most the same fit quality as is shown in Fig. 01 Simulation 
of the Stokes I and V profiles in the regions of these 
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lines results in almost the same configuration of the mag- 
netic field structure as obtained from the analysis of the 
stronger Fe n lines, for which all four Stokes parameters 
were taken into account during the simulation. 



o 




■1500 ' — ' — ' — ' — ' — 1 — ' — ' — ' — ' — 1 — ' — ' — ' — ' — 1 — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — 1 
-0.25 0.25 0.5 0.75 1 1.25 

Phase 

Fig. 6. Longitudinal magnetic field data with phase ob- 
tained from analysis of Fen A6432.68 (solid line), Fen 
A5018.44 (dashed line), Fen A6516.08 (dotted line) and 
Fen A5197.58 (dash-dotted line). The open squares with 
la error bars correspond to the LSD data published by 
Wade et al. (20QQb|). 



4.3. Integral magnetic field characteristics 

In order to check the agreement of the derived magnetic 
field model with other available magnetic field data for 78 
Vir, we calculate the intensity-weighted, averaged (over 
the visible stellar disk) longitudinal magnetic field B\ and 
the normalized equivalent widths of the Stokes Q and U 
profiles for all the analysed phases. The normalization pro- 
cedure is performed in accordance with the method de- 
scribed by Wade et al. il2t)l)l)bi over the passband of the 
analysed line profiles. 

All the available longitudinal magnetic field mea- 
surements for 78 Vir are plotted in Fig. of 
Leone & Catanzaro H2001[l . This figure shows the 
good agreement of the B\ measurements obtained 
by Wade et al. (2000b) using the Least-Squares 
Deconvolution (LSD) technique (Donati et al.[l997) with 
the majority of other observational data. Therefore, in this 
paper we just compare our results with the LSD longitu- 
dinal field data. 

As demonstrated by Wade et al. (|2000b|) . the 78 Vir 
Stokes Q and U LSD equivalent widths are proportional 
to the BBLP measured at that phase (LeroyEESI) with a 
line scaling factor. This factor might be different for the 
analysed lines (see the last column in Table |2J) , but is the 
same for the both Stokes Q and U equivalent widths for 
a particular spectral line. 

Neither the B\ measurements nor the BBLP data are 
included directly into the model minimization procedure. 
However, given that the model successfully reproduces the 
Stokes profiles with which these quantities are fundamen- 
tally related, we should expect that the model magnetic 



0.1 




Phase 



Fig. 7. Comparison of scaled net linear polarization data 
obtained from the best fit simulations of Fen A4923.927 
(solid line) and Fe II A6516 (dashed line), and from LSD 
profiles (Wade et al.|2000b| (open squares) for 78 Vir with 
respective la error bars. 

field configuration is capable of reproducing them. Using 
the same stellar atmosphere model, the longitudinal field 
calculated from the model fits to the strongest Fen line 
A5018.44A and to the weak Fen line A6432.68A show a 
good agreement with the LSD B\ data (see Fig. EJ. The 
longitudinal fields calculated from the simulation results 
for the other Fe n lines appear to be shifted downward in 
longitudinal field intensity relative to the LSD Bi variation 
(see Fig. EJ). 

The calculated Stokes Q and U equivalent widths for 
all the analysed line profiles are multiplied by the respec- 
tive line scaling factor in order to fit them to the LSD 
Stokes Q and U data. For each analysed line a scaling 
factor is determined with the help of the Least Squares 
method. The respective results are given in Tabic |2J 
taking into account that the LSD equivalent widths of 
Stokes Q and U profiles have a line scaling factor of 0.1 
(Wade et al. l2000b|l from comparison with the BBLP data 
CLerov Hlj95|l . Fig. [7] presents the Stokes Q and U equiv- 
alent widths variability with phase, derived from the best 
fit simulations of Fen A4923.927 (with Q = 109°) and Fen 
A6516.08 (with fl = 128°). They are scaled by a factor of 
0.08 (for A4923) and 0.03 (for A6516) and plotted over 
the Stokes Q and U equivalent widths derived from LSD 
profiles. It is remarkable that the simulated polarimetric 
data for Fen A4923.927A and A5018.44A lines have al- 
most the same scaling factor as the LSD Stokes Q and 
U equivalent widths. Theoretically, the line scaling fac- 
tor depends on the Lande factor and on the degree of line 
saturation (Wade et al. l2000a|l . Tabled shows that in gen- 
eral the lines with a high mean Lande factor have also a 
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Table 4. The same as at Table |3 but for Cm and Tin 

lines 



Line 


Or II 


TT 

Or II 


TT 

Or II 


1 .. T T 

Or II 


t~*~ TT 

Or II 


tt 

Or II 




Ti II 


Ti II 




A 


4592 


4634 


5237 


5310 


5407 


5421 


O'er 


5188 


5336 


O'er 


2 

Xi 


14.82 


13.12 


10.51 


6.47 


6.20 


9.48 




21.53 


12.45 




A 2 


8.83 


6.39 


12.70 


8.64 


11.11 


11.50 




24.28 


9.30 




ft 2 
ft 2 

6Xu 


5.62 


6.18 




8.18 


8.55 


11.76 




8.77 


5.79 




6.48 


5.88 




6.63 


9.17 


12.33 




7.45 


6.45 




v 2 
Aw 


Q OK 

o.yo 


1 on 
1 .89 


1 1 ftn 
11. OU 


*7 a n 
I .4y 


O Tft 
8.(0 


11. Z ( 




1 r F1 

10. ol 


o £n 
S.oU 




,,2 
X 


A Ta 
4. lb 


4.1s 


ft O A 

0.84 


I. It 


o no 
2.98 


a nn 
4.09 




T c o 
f .58 


a on 
4.20 




CJr, kG 


175 


186 


185 


182 


161 


171 


30 


287 


271 


40 


oi, 10 -3 


10.2 


9.5 


9.7 


8.9 


9.6 


8.6 


1.5 


3.7 


3.8 


1.4 


Ai 


19° 


1 r 

15 


16 


25° 


17° 


21° 


r o 

5 


28° 


31° 


17° 


Si 


-41 


-42° 


-45 


-49° 


-40° 


-46° 


r o 

5 


-39° 


A 1 O 

-41 


19° 


02, 10 -3 


4.4 


4.1 


3.8 


3.9 


3.6 


3.7 


0.3 


2.8 


3.5 


1.0 


>2 


103° 


120° 


111° 


115° 


109° 


120° 


no 

8 


113° 


113° 


19° 


5 2 


-17° 


-18° 


-17° 


-12° 


no 

-2 


-19° 


10° 


oo 

-8 


1 AO 

-14 


i no 

12 


Q 


109° 


117° 




123° 


115° 


109° 


14° 


97° 


90° 


19° 




26° 


26° 


26° 


26° 


22° 


27° 


5° 


25° 


18° 


5° 


log(El/N tot ) 


-3.90 


-4.02 


-4.04 


-3.91 


-4.18 


-4.17 


0.19 


-5.48 


-5.68 


0.22 


Vr 


-8.25 


-8.02 


-8.80 


-8.41 


-8.05 


-8.25 


1.23 


-5.61 


-7.91 


1.26 


V e sin i 


11.7 


10.8 


11.8 


11.5 


10.9 


12.2 


1.1 


13.4 


13.3 


1.5 


P 


123° 


119° 


126° 


131° 


127° 


124° 


5° 


117° 


112° 


5° 


a ,10- 3 


6.1 


5.2 


5.5 


5.2 


5.1 


4.9 


1.4 


2.5 


2.9 


1.7 


a,10 -3 


5.0 


5.1 


5.0 


4.6 


5.1 


4.4 


1.4 


2.1 


2.2 


1.7 


Bp, kG 


3.52 


3.86 


3.71 


3.36 


3.32 


3.05 


0.4 


2.46 


2.41 


0.5 




-11° 


-10° 


-11° 


-8° 


-8° 


-7° 


8° 


-18° 


-23° 


8° 


5 P 


-33° 


-29° 


-36° 


-41° 


-37° 


-34° 


10° 


-27° 


-22° 


11° 


B», kG 


-3.44 


-3.78 


-3.62 


-3.29 


-3.24 


-2.99 


0.4 


-2.45 


-2.41 


0.5 


A„ 


169° 


170° 


169° 


172° 


171° 


173° 


8° 


161° 


156° 


8° 


Sn 


33° 


29° 


35° 


40° 


37° 


34° 


10° 


27° 


22° 


11° 



comparatively high line scaling factor. Nevertheless, this 
dependence is not clearly pronounced by the sample of 
derived data. This is because, for some lines, this effect 
is reduced by blend contamination of the analysed Stokes 
Q and U profiles (for example the Fen A5100 line group, 
A5197.03 and A5362.87 lines) or by weak Stokes Q and 
U profile variability, which results in a comparatively low 
precision of the calculated equivalent widths. 

The LSD net linear polarisation (in the same way 
as the BBLP) reflect the contributions of many spectral 
lines with significant polarization. These lines belong to 
a number of chemical elements, which apparently have 
non-uniform and usually quite different abundance distri- 
butions. Such a composition of polarised features results 
in LSD (or BBLP) data that may differ substantially from 
the Stokes Q and U equivalent widths derived for a par- 
ticular Fen line (see Fig. EJ). This is illustrated by the 
rather large differences between the two simulated curves 
in Fig. Moreover, lines of elements with non-uniform 
abundance distributions can provide substantially differ- 
ent values of the sky-projected position angle of the stel- 
lar rotation axis (see Sect. I4.4|l . and can exhibit totally 
different net linear polarisation variations (e.g. 180° out 
of phase with those shown here). The spectral resolution 



and the observational errors of the available spectra do 
not allow us to map confidently the detailed distribution 
of chemical abundances nor the local patterns of the sur- 
face magnetic field. The uncertainties in the details of the 
abundance distributions and the magnetic field structure 
limit the precision of f2 as well as the quality of simulated 
Stokes Q and U equivalent widths. Given these limita- 
tions, we would characterise the agreement of the observed 
and simulated variations shown in Fig. [7] as very accept- 
able. 

4.4. Lines of other elements 

78 Vir is a chromium-rich star (Cowley et al. I1969f) 
and has some strong Cm lines, which are located 
in the observed spectral range. The most promi- 
nent of them, Cm A4592.049, A4634.07, A5237.329, 
A5310.687, A5407.604 and A5420.922A, have been se- 
lected for simulation. Atomic parameters of spec- 
tral lines were extracted from the VALD database 
(Kupka et al. 1999) and from Raassen & Uylings 1|1998[) 
(ftp://ftp.wins.uva.nl/pub/orth) in the case of Cm lines. 
These lines are not significantly blended by lines of other 
elements and show variability of the Stokes IV param- 
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c) Q/h 



d) U/h 
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Fig. 8. The same as at the Fig. but for Cm A5310.687A. Non-uniform chromium distribution in the stellar atmo- 
sphere is assumed during this simulation (dashed line). For comparison we show also the best fit simulation result 
(\ 2 = 4.24) for uniform chromium distribution (dotted line). 



eters with phase. Some show marginally-detected Stokes 
QU signatures as well. 

We analyse also the two strong Tin A5188.68 and 
A5336.771 lines. These Ti lines are not significantly 
blended (although we take into account a small contri- 
bution of the Vi A5188.885 blend to the profile formed by 
Tin A5188.68). 

Table 0] presents the results of the best fit simulations 
for the Cr and Ti lines. The best fit quality is on average 
similar to that obtained for the Fen lines. The Stokes I 
profiles for these particular lines vary weakly with rota- 
tional phase, suggesting the presence of non-uniform dis- 
tributions of Cm and Tin in the atmosphere of 78 Vir. 
To take into account the non-uniform abundance distri- 
butions the stellar surface is divided into 30 areas, each 
characterised by an independent local abundance of the 
analysed chemical element. The local abundances are in- 
cluded as free parameters in the simulation procedure and 
the model operates in this case with 41 parameters. The 
best fit simulations of the aforementioned lines tentatively 
suggest that chromium is enhanced in the vicinity of the 



negative magnetic pole, while titanium is underabundant 
in this region (see Fig. 110(1 . The improvement in the fit 
to Cr II A5310 resulting from a non-uniform abundance 
distribution is highly significant (reduction of xj by over 
40%), although the improvement for Ti n A5336 is much 
less so (just 10%). The averaged chromium and titanium 
abundances are given in Table 

It is remarkable that the best fit simulations of Cm 
lines provide almost the same global magnetic field struc- 
ture and the sky-projected position angle of the stellar 
rotation axis as the Fe lines (see Tables 0J) . Meanwhile, 
the Ti II simulation results in somewhat lower values of Q 
and (3. These differences may results from the non-uniform 
Ti II distribution in the stellar atmosphere. 

5. Discussion 

The abundance of Fen derived in this study is 
\ogFe/N tot = -3.16 ±0.20 (see Table EJ. From analy- 
sis of spectra of 78 Vir in the regions 3850-3870 A and 
3868-4650 A Adelman Ijl973b|l obtained an Fen abun- 
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Fig. 9. The same as at the Fig.|2 but for Ti n A5336.77lA. Non-uniform titanium distribution in the stellar atmosphere 
is assumed during this simulation (dashed line). For comparison we show also the best fit simulation result (x 2 = 4.58) 
for a uniform titanium distribution (dotted line). 



dance of log Fe/Nn = — 2.89dex. The differences in the 
iron abundances may be due to the inclusion/exclusion of 
polarized transfer, microturbulence, effects of Balmcr line 
wings, etc., as well as from the different effective tempera- 
tures that have been employed for the abundance analysis 
(for 78 Vir Adelman l|1973a|l applied T cS = 9950K, sig- 
nificantly higher than the temperatures employed in our 
model. 

For the 11 different Fe features studied here, we ob- 
serve no systematic difference in strength between weak 
and strong lines. Therefore, the results of this study are 
consistent with the absence of important Fe stratification 
in the atmosphere of 78 Vir. 

According to our results the global magnetic field 
structure of 78 Vir is well-described by a slightly decen- 
tered magnetic dipole. Usually the starting point of a sim- 
ulation (in the free parameter hyperspace) is chosen to cor- 
respond to a central magnetic dipole, but for all the anal- 
ysed lines the final model results in a slightly decentered 
magnetic dipole. This fact is in good qualitative agree- 
ment with the results previously obtained by Borra (Q980 ). 



In that paper Borra considered the classical decentered 
dipole model, where the magnetic dipole size is insignifi- 
cant in comparison with the stellar radius, and obtained 
ao—0.2. Our model provides a smaller a = 0.006 ± 0.002 
due to the non-zero dipole size (2a = 0.008 ± 0.003), that 
allows us to take into account the non-symmetrical (with 
respect to the stellar center) magnetic field configuration 
of the star. The magnetic field intensity and location of 
the positive and negative magnetic poles in the stellar ro- 
tational reference frame are: 



R, = 3.4±0.7 kG, X p 



-9°±5 C 



-33°± 12° 



B„ = -3.3±0.7 kG, A n = 171°±5°, <5„ = 33°± 12°. 



(7) 



Due to the shift of the dipole from the stellar centre, the 
positive magnetic pole has a slightly stronger field inten- 
sity than the negative pole, but not all the analysed lines 
are equally sensitive to this difference. Nevertheless, the fi- 
nal model is hardly distinguishable from a symmetric cen- 
tral dipole model, which provides a slightly poorer agree- 
ment between the simulated and observed data. For Fe n 
lines the centered dipole model results in a x 2 " mnc tion 
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Fig. 10. Illustration of the a) chromium and b) titanium abundance distribution in the stellar atmosphere of 78 Vir. 
Location of the positive rotation pole is specified by open circle. The stellar image has been rotated in the plane of 
the page by the positional angle 0=110°. 



which exceeds by 3%— 5% the best fit result from the de- 
centered dipole model. This fact can be related to the 
quality of the data and to the model sensitivity. The dipole 
offset could be a result of the particular geometry of 78 
Vir, because we never completely see the positive mag- 
netic pole of the star. The modeled global magnetic field 
configuration of 78 Vir is illustrated in Fig. ^2 according 
to the best fit results obtained for Fell A5018A. Five of 
the Fen lines analysed using the entire Stokes vector re- 
sult in a plane-of-sky orientation of the rotational axis of 
Q = 110° ± 17°, while the last two lines (Fe n A6432A and 
A6516A), with comparatively weaker variability of the ob- 
served Stokes Q and U profiles, provide higher values of 
this angle. Fig. ^| shows the dependence of the Xq an d 
Xu functions on the angle for the three strong Fe n lines 
(see Subsec. r4.2.1|) . Both functions reach their minima at 
SI = 109° and SI = 289°. The second value describes a 
model which is indistinguishable from the configuration 
specified in Table [3] due to the definition of the Stokes Q 
and U parameters. 

As discussed in Sect. 3, 78 Vir is a probable member 
of the Ursa Major stream. It approaches us with a mean 
radial velocity V T = -8.1±1.0 km s" 1 (Wade et al. 2000a). 
This velocity varies by about ±2 km s _1 with the phase 
of stellar rotation ('Preston ll969l this work). The possible 
variability of the radial velocity is not taken into account 
in the simulation procedure. Respectively, our estimations 
of the radial velocity for different Fe II lines are distributed 
around the aforementioned value within the range ±1.5 
km s _1 . 

Nevertheless, the V r variability through the rota- 
tional cycle has been determined from the strong 
Fen lines A4923.927A A5018.44A and A5169.033A (see 



Subsec. 14.2.1(1 and for the other Fell lines as well. We 
show that the mean radial velocity derived from the simu- 
lation varies from line to line (see Table|2J). The other Fen 
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Fig. 11. Illustration of the surface magnetic field geome- 
try of 78 Vir. The location of the positive rotation pole is 
indicated by the open circle. We always see the negative 
magnetic pole and only for phase (p = 0.0 are the positive 
and negative magnetic poles visible simultaneously. The 
stellar image has been rotated in the plane of the page by 
the positional angle Sl=110°. 



lines provide a similar A V r variability. The precision of the 
mean radial velocity determination depends on the obser- 
vational errors and on the quality of the description of the 
surface magnetic field structure. For the blended profiles 
it also depends on the accuracy of oscillator strengths. 

The analysed lines show weak, coherent variations of 
AV r with the period of stellar rotation (see Fig. 0J. This 
suggests a possible moderately non-uniform distribution of 
Fell. Nevertheless, the non-uniform iron distribution will 
not lead to a significantly different surface magnetic field 
structure. The simulation of Cm lines with the assump- 
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Fig. 12. Dependence of a) Xq anc ^ b) X\j on the angle calculated for three strong lines Fen A4923.927A (solid 
line), A5018.44A (dashed line) and A5169.03A (dotted line). The calculations were performed under the condition of 
stability of all the best fit parameters for these lines (see Table except the fi. 



tion of a non-uniform chromium distribution resulted in 
almost the same field structure that we obtained from the 
simulation of Fen lines. 

The other best fit parameters arc the rotational axis 
inclination i = 24° ± 5° and the dipole axis obliquity 
/3 — 124° ± 5°, which are in a good agreement with the 
estimates of Leroy et al. i[199fi|l . The derived 14 sin i = 
12 ± 1 km s^ 1 provides V c — 29 ± 4 km s" 1 and is consis- 
tent with that of Preston l|1971|l . 

The derived surface magnetic field variability inter- 
val ranges from 2.1 kG to 3.2 kG and covers the value 
B s =2.9 kG estimated by Preston l|1971|) . Besides, as Fig.© 
shows, the longitudinal magnetic field variation obtained 
from the two Fell lines simulation is also in good agree- 
ment with the LSD B\ data (Wade et al. 12000 bjl . These 
facts justify the applicability of the slightly decentered 
magnetic dipole model with the aforementioned values of 
free parameters for the global magnetic field structure de- 
scription at 78 Vir. 

Unfortunately in the case of 78 Vir the intensity of 
the Stokes Q and U profiles is similar to the observational 
errors. The data are unsuitable for performing an analysis 
of the small-scale structure of the magnetic field, using a 
more sophisticated technique such as Magnetic Doppler 
Imaging (MDI; Kochukhov & Piskunov 12002(1 . However, 
it appears, given the relatively good agreement between 
the intensity, structure and variability of the observed and 
simulated Stokes QU profiles, that the real magnetic field 
of 78 Vir does not depart strongly from the configuration 
derived here. At the same time, we have noted that the 
Stokes QU profiles are not, at some phases, fit to within 
their errors. Therefore these data, notwithstanding their 
relatively low S/N, already suggest the limitations of the 
MCD model framework. 

A similar analysis of the Stokes IV profile variability 
was performed by Kochukhov et al. {5002) for a 2 CVn. 



Those authors applied MDI using "multipolar regulariza- 
tion" , analysing the Stokes IV line profiles and obtained a 
good agreement between observed and computed profiles 
for Fen, Cm and Sin lines. The quality of the Stokes I 
and V profile fits obtained in this paper is only marginally 
poorer than that obtained by Kochukhov et al. (2002), al- 
though remarkably the magnetic field model framework 
employed here is much simpler. 

New spectro-polarimetric observations of 78 Vir with 
higher spectral resolution and signal-to-noise ratio (espe- 
cially for the Stokes Q and U spectra) are required for 
further refinement of the global magnetic field structure. 
With such data, a model with additional parameters (such 
as MDI) could be applied, in order to study the local mag- 
netic field topology and the detailed relationship between 
the magnetic field and the abundance distributions. 
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